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ABSTRACT 

The results of a search for sub-km Kuiper Beh Objects (KBOs) with the method of serendipitous 
stellar occultations are reported. Photometric time series were obtained on the 1.8m telescope at the 
Dominion Astrophysical Observatory (DAO) in Victoria, British Columbia, and were analyzed for the 
presence of occultation events. Observations were performed at 40 Hz and included a total of 5.0 
star-hours for target stars in the ecliptic open cluster M35 (/3 = 0.9°), and 2.1 star-hours for control 
stars in the off-ecliptic open cluster M34 (/3 = 25.7°). To evaluate the recovery fraction of the analysis 
method, and thereby determine the limiting detectable size, artificial occultation events were added to 
simulated time series (1/f scintillation- like power-spectra), and to the real data. No viable candidate 
occultation events were detected. This limits the cumulative surface density of KBOs to 3.5 x 10^° 
deg^^ (95% confidence) for KBOs brighter than mfl=35.3 (larger than ~860m in diameter, assuming 
a geometric albedo of 0.04 and a distance of 40 AU). An evaluation of TNO occultations reported in 
the literature suggests that they are unlikely to be genuine, and an overall 95%-confidence upper limit 
on the surface density of 2.8 x 10^ deg~^ is obtained for KBOs brighter than m/j=35 (larger than 
~1 km in diameter, assuming a geometric albedo of 0.04 and a distance of 40 AU) when all existing 
surveys are combined. 

Subject headings: occultations, Kuiper Belt, Solar System: formation 



1. INTRODUCTION 

The physical and orbital properties of the Kuiper Belt 
Objects (KBOs) are believed to provide valuable infor- 
mation about the formation of the solar system. To date, 
directly observed KBOs include objects with dia meters 
rangi ng from ^25km (eg. 1999 DAg, [C ladma n et al.l 
I2OOII: and 2003 BH91, Bernstein et"alll20b4f) to 2400km 
(eg. Eris, [Brown et al. 2006 ). KBOs with diameters be- 
low a few tens of kilometers are beyond the detection of 
current ground and space-based telescopes. 



The KBO size distribution is the outcome of accre- 
tion and coll i sional evolution in the o u ter sola r system 



(IStern 


119961: iDavis & Farinellal 119971 IStern & ColweUl 


19971: 


Kenvon & Luul I1999U iKenvon & BromlevI 12004 


Pan & Sarill2005l). The cumulative luminositv function 



(CLE) of the KBOs provides an indirect means of deter- 
mining their size distribution, and is observed to have 
the form: N{< ttir] tur < 25) = I0"(™«-^o). The 
KBO size distribution is then a power-law of the form: 
dA'^ oc r~''dr, with slope —q, where q is related to a via: 
g = 5a + 1 (Gladman ct al. 2O 0..ll, GOl)- Cur rent obser- 
vations indicate a ss 0.68 (jFraser et al.ll2007[). or equiv- 
alent ly , q = 4.6. HST observations bv iBernstein et ahl 
(|2004t ) suggest a lower CLE slope for niR > 25. 
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Current theoretical models of the accretion and frag- 
mentation processes in the early solar system suggest 
that q remains constant only for larger objects, and 
that the slope of the size dist ribution changes at t he 
'break radius', rf, <100 km (iKenvon fc Luul Il999al lbl: 
IKenvon fc BromlevI [200I iPan fc Saril 120051) . Objects 
smaller than the break radius would be in collisional 
equilibrium and their size distribution would h ave the so- 
called Dohnanyi equilibrium slope of g « 3.5 (jDohnanvil 
Il969| ). Object densities and bulk-moduli can influence 
the slope for the small objects, and th e position of the 
break-radius (jKenvon fc Bromlevll2004D . 

The deepest dir ect observations to date 
([Bernstein et al.l |2004J ) were performed with HST 
and reached a limiting magnitude of R=28.5 (corre- 
sponding to an object diameter of '^20km for an albedo 
of 0.04 at 40AU). The number of observed objects 
was ~25x lower than the extrapolation of the current 
CLE slope. This absence of faint objects indicates a 
decrease in the density of the Kuiper Belt beyond 50 
AU; and a dearth of sub-lOOkm KBOs, in general. 
Direct observation is not a feasible option for pursuing 
the much fainter km-sized KBOs. 

Some indirect methods of inferring the slope in the 
range of the km- s ized objects have been explored. 
IStern fc McKinnonI ()2000D used Voya ger 2 images of Nep- 
tune's moon Triton to analyze cratering and infer a slope 
of g « 3 for the size distribu t ion o f the km-sized im- 
pactors. IKenvon fc WindhorstI (|200lD have placed upper 
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limits on small object densities by examining the optical 
and FIR sky surface brightnesses of the ecliptic. They 
found that a straight extrapolation of the CLF to R«40- 
50 mag would yield an ecliptic surface brightness higher 
than that observed. 

One promising method of indirect detection for the 
smaller KBOs is tha t of ser e ndipit o us stellar occulta- 
tions (iBailevI [1976 1: iDysonI fl99^ 1 rown fc Websted 



1997t iRoaues fc M oncuguet 20 00l: iCooravl 120031: 



Coorav &: Farmed 120031: " iGaudil ~ l2004t iNihei et all 
200711 ■ A km-sized KBO passing through the line of 



sight to a star will partially obscure the light from that 
star to reveal the presence of the otherwise invisible 
object. Such an occultation cannot be predicted and the 
practice relies upon serendipity. 

A 1 km diameter KBO at 40 AU subtends and angle of 
only 34 /uas, and the lightcurve resulting from an occul- 
tation by such an object is expected to exh ibit a series of 
incre ases and decreases due to diffraction (jRoques et al.l 
Il987t ) (henceforth R87). At opposition, the orbital veloc- 
ities of the Earth (30 km s^"'^ ) and a typical KBO at 40 
AU (^4 km s^^ ) produce a retrograde relative velocity 
of V w26 km s^ for the occulter, and the characteristic 
diffraction shadow cast by a 1 km KBO will pass over a 
terrestrial observer in a fraction of a second. Observa- 
tions must be performed at >40 Hz to sample the flux 
oscillation in an observed lightcurve. 

We have performed a KBO occultation search with the 
1.8m Plaskett Telescope at the Dominion Astrophysical 
Observatory (DAO) using a 40Hz camera. Models of oc- 
cultation lightcurves are presented in Section [2l Obser- 
vation details are given in Section [31 with data reduction 
and analysis methods described in Sections [4l and [5l re- 
spectively. Our results are presented in Section [H and 
reported KBO occultation detections are evaluated in 
Section [7] A surface density upper limit is presented in 
Section [8l with a discussion and summary provided in 
Sections [H and [TOl respectively. 

2. MODELING OCCULTATIONS 

The object diameter at which diffraction effects begin 
to dominate a KBO occultation lightcurve is known as 
the Fresnel scale, p = \z/2, where A is the wavelength 
and z is the distance between the observer and the occult- 
ing body. For a KBO occultation (2;=40AU) observed in 
visible light (550nm), the Fresnel scale is 1285m. Fresnel- 
Kirchhoff diffraction theory should be used to model oc- 
cultation lightcurves for objects with sizes of order 1 Fres- 
nel scale unit (Fsu). Following R87, we have modeled oc- 
culters as circular disks, with the masking disk assembled 
by placing progressively smaller rectangular masks at the 
periphery of the structure. Five 'orders' of such rectan- 
gles were used for our models. An analytic circularly- 
symmetric solution to the Fresnel-Kirchhoff equation ex- 
ists (R87), but we chose to develop our KBO disks from 
assemblies of rectangles to build in flexibility of object 
shape for future work. 

2.1. Passbands and Source Dimensions 

Intensities were numerically integrated over a range of 
wavelengths to model a broadband observation of an oc- 
cultation. This was accomplished by averaging evenly- 
weighted diffraction patterns at wavelengths distributed 



uniformly across the passband. Lightcurves for 550nm 
monochromatic light and 400-700nm broadband light are 
compared in Figure [1] The extended 'ringing' present in 
the monochromatic pattern is attenuated by interference 
in the broadband pattern; only the first peak away from 
the pattern center is well preserved. 
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Fig. 1. — Diffraction profiles for monochromatic light {550nm) 
and broadband visible light (400-700nm). The occulting object 
was 250m in radius at a distance of 40 AU. 

It was unnecessary to weight the intensities to match 
the spectral energy distribution (SED) of the target star, 
the spectral throughput of the telescope, and the spectral 
sensitivity of the CCD. Even-weighting offered greater 
precision than needed for work with photometric time 
series having > 1% RMS variability. 

The target stars were modeled by averaging evenly- 
spaced point sources over the projected stellar disks. A 
point source density of 100 Fsu~^ was used (~100 points 
for the stellar disks of the M35 targets projected at 40 
AU) . This density was more than sufficient to provide an 
accurate model of the shadow , consi stent with models 
found in IRoaues fc Moncuguet (|2000f ). 

2.2. Conversion of the Shadow Profile to a Time Series 

Profile 

Intensities for the diffraction-dominated occultation 
shadows were computed as a function of position with 
respect to the center of the projected shadow. The 
relative velocity was used to convert positional coordi- 
nates to time coordinates to determine the corresponding 
lightcurve that would be recorded by an observer passing 
through such a shadow. 

An object observed at a given ecliptic latitude must 
have an orbital inclination equal to, or greater than, 
that latitude. As the fraction of TNOs is observed to 
decrease with incre asing orbital inclination (jJewitt et al.1 
119961 : [Brownll200lT ): an observed object is most probably 
very near its highest latitude, where its velocity vector 
lies nearly parallel to the ecliptic plane. As only the 
ecliptic latitude of observation is known, the orbital in- 
clination was taken to be equal to the ecliptic latitude. 
This is the most probable configuration, and it is illus- 
trated in Figure [3 

The relevant velocity is the component of the relative 
velocity projected onto the plane perpendicular to the 
line-of-sight. Letting 9 represent the angle between the 
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Fig. 2. — The orbital geometry used to derive the relative ve- 
locity of a body orbiting on an inclined orbit. Vectors labeled 'r' 
and 'v' denote heliocentric distance and velocity, respectively. The 
vector 'd' is the distance from the Earth to the occulting KBO. 
The subscript 'k' denotes parameters referring to a KBO, and 'e' 
denotes those referring to the Earth. The subscript 'p' is used 
to indicate quantities projected onto the ecliptic. The Earth and 
KBO are indicated with ©, and • symbols, respectively. The an- 
gle £ is the solar elongation projected onto the ecliptic plane. The 
angles a and /3 are measured in the ecliptic plane; ol measures the 
angle between the position vector of the Earth, re, and the line-of- 
sight to the target star projected onto the ecliptic plane, r,p (r*p 
is parallel to dp). 

line-of-sight, d, and the relative velocity vector, Vrel^ that 
component is Vrei-perp = \vrei\ sin{9) . Calculation of 9 
can be performed with the dot product d ■ v^ei ■ 

9 = arccos i^TTp^A ■ (1) 

The vectors d and v^ei were expressed in Cartesian coor- 
dinates in terms of the quantities displayed in Figure [2l 
The quantities: Vkp (heliocentric distance to the KBO 
projected onto the ecliptic place); angles /3, and a; and 
dp (observer-to-KBO distance projected onto the eclip- 
tic) were determined. 

|rfep| = |rfc| cos(i), (2) 

Of = TT — £, (4) 

|dp| = + keP - 2 \vkp\ |re| cos(a - /3). (5) 

The position coordinates are: 

(xe, ye, Ze) ^ (|re|, 0, 0), and (6) 
/|rfep| cos(a - /3)\ 

{xk, Vk, Zk) = |rfep|sin(a - /3) . (7) 
V |rfe|sin(i) / 

The velocity components are: 

[Vxe, Vye, V^e) = (0, |Ve|, 0), and (8) 



/-|vfc|sin(«-/?)\ 
{vxk, Vyk, Vzk) = [ |vfe| cos(a - /3) . (9) 
V / 

The vectors d, and Vrei are then: 

d = {xk-Xe, Vk-TJe, Zk-Ze), and (10) 

Vrei = (Vxk - Vxe, Vyk - Vye, V^k - V^e) (H) 

Vectors d and Vrei were substituted into Equation [T] to 
compute 9] and, in turn, the perpendicular component 
of the relative velocity Vrei-perp = Wrei \ sui{9). 

If r/cp is less than re, the algorithm will fail for cer- 
tain elongations as they no longer exist. For example, 
an object on an orbit with r^p < cannot be observed 
at opposition; its projected position is in the same hemi- 
sphere as the sun. 

3. OBSERVATIONS 

The open clusters: M35 (/3 = 0.9° N), and M34 
(/3 = 25.7° N) were chosen as our primary target field, 
and our off-ecliptic control field, respectively. Two adja- 
cent stars were monitored in each field. The properties 
of the targets are listed in Table [TJ Values of my, B- 
V, and the spectral (MK) typ es for the M34 stars were 
taken from iKharchenko et a L (2004'). Those for M35 
were taken from ISung fc Bessell (:,199a) , but the spec- 
tral types were not available and had to be inferred by 
comparison of My and B-V values to those for other 
upper ma in sequence stars of kno wn spectral type (re- 
ported in IKharchenko et all I2004D in the cluster. Val- 
ues for the stellar r adii, R*, w ere estimated based on 
their spectral types (|Co:>3 12000( ). The distance moduli 
(|Saraiedini et al.l |2004[ ) were used to compute the dis- 
tances with (m — M)y = 5 log(c?) — 5. Values of i?40Au 
denote the stellar radii projected at 40 AU. Both clus- 
ters are sufficiently distant that the selected stars have 
projected diameters less than the Fresnel scale for visible 
Hght at 40AU (~1300m). 

Observations were conducted in 2-3"seeing conditions 
with the 1.8m Plaskett Telescope at the Dominion As- 
trophysical Observatory (DAO) on the UT nights of 
2004 January 2, and 2004 January 3. These dates were 
chosen to place the primary target near an elongation 
of £ = 180°, where occulta tion probability is highest 
(jRoques &: Moncuquetl 12000^. On those dates, the so- 
lar elongations of M35 and M34 were e = 169°, and 
e = 128°, respectively. 

A Texas Instruments TC253-SPD CCD was operated 
with a San Diego State University 2 (SDSU2) controller, 
and data were binned 8x8, and sub-rastered down to a 
smaller region of the chip to enable the CCD to be read- 
out at 40IIz. The binned pixel scale was 2.67" pixel"^. 
Images taken early in the observing run were 30x24 pix- 
els, and images taken later were 25x16 (8x8 binned, 
unprocessed sizes). One hundred images were taken in 
repeated sequences, and were stored in 3D-FITS format. 
Each was followed by a ^0.5s overhead as the controller 
initiated the next 100-frame sequence. In total, 357,000 
images were obtained for M34, and 794,000 were ob- 
tained for M35. 

4. DATA REDUCTION 

A bias frame was constructed by median-combining 
3600 25ms dark exposures. The count level was ~3000 
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TABLE 1 

Properties of the stars in the target field (M35), and 

CONTROL FIELD (M34). 

M34 Targets" 



Position 


East (left) 


West (right) 


ID 


TYC-2853-22-1 


TYC-2853-679-1 


RAj2000 


02:41:58.31 


02:41:56.78 


Decj2ooo 


+42:47:29.7 


+42:47:22.4 


/3 


25.7° N 


25.7° N 


my 


8.360 


8.444 


B-V 


0.002 


-0.005 


MK 


B9V 


B9Vp Mn 


R* 


3.0 Rq 


3.0 Rq 


-'^40AU 


640m 


640m 


(m-M)y 


8.98±0.06 


8.98+0.06 


Dist. 


625±17 pc 


625+17 pc 


M35 Targets 


Position 


East (left) 


West (right) 


ID 


TYC-1877-704-2 


TYC-1877-266-1 


RAj2000 


06:09:20.75 


06:09:18.50 


Decj2ooo 


+24:22:08:6 


+24:22:05.4 


/9 


0.9° N 


0.9° N 


my 


10.656 


9.764 


B-V 


0.178 


0.167 


MK 


B9V 


B9V 


R* 


3.0 Rq 


3.0 R0 


-'^40AU 


360m 


360m 


(m-M)y 


10.21+0.12 


10.21+0.12 


Dist. 


1101+60 pc 


1101+60 pc 



ADU with an estimated sky contribution of ~2 ADU. 
One thousand 100ms sky-flat cahbration images were 
median-combined to produced a flat which showed no 
structure and had pixel-to-pixel variabiUty of less than 
1%. The true variability across the field was, therefore, 
assumed to be negligible, and no flat-field correction was 
applied. Examples of raw and pre-processed images are 
shown for the M35 field in Figure [31 
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Fig. 3. — Unprocessed (left) and a pre-processed (right) M35 
frames. The frames are 25x16 pixels, and 21x12 pixels, respec- 
tively. With 8x8 binning, the processed image is 56" X 32" . 

As the high-speed camera was installed in place of the 
guider, guiding was unavailable during the observations. 
Without guiding, the Plaskett telescope's tracking mo- 
tor introduced an oscillation in right ascension with an 
amplitude of A fa 8 arcsec, and a period of T fnA min- 
utes. This corresponds to a maximum drift speed across 
the chip of 2ttA/T «0.2 arcsec s~^, or 0.08 pixels s~^. 
The target star's motion between pixels was too small to 
induce a calibration error on the sub-second time scale 
of a KBO occultation. 

Photometry was per formed with SExtractor v2.3 
(|Bertin fc Arnoutsl[l996l ) . The stars in the field were not 
crowded, and aperture photometry was performed with 
5 pixel (13.3 arcsec) apertures (diameter). Other SEx- 
tractor parameters were set at their default values. Flux 
values for the targets in the 100-frame sequences were 



concatenated into longer sequences of ~50,000 - 100,000 
frames. The 0.5s gaps occurring between the sequences 
were masked and ignored. The gaps are cosmetically un- 
appealing, but only diminish the usefulness of the data 
in proportion to the duty cycle reduction they represent 
(- 17%). 

Each time series was normalized to have a constant 
mean flux. A sufficiently large number of simultaneously 
observed targets would make it possible to normalize by 
dividing by the summed flux of the other targets. With 
two targets of comparable brightness, this method would 
have increased the noise level in a time series by ~40% 
(•\/2). Each time series was therefore normalized by di- 
viding by a Gaussian-smoothed version of itself. To con- 
struct smoothed time series, each point was replaced by 
the Gaussian-weighted average of the flux measurements 
preceding it and following it, such that for flux /, the 
smoothed flux at the j*'' time coordinate was: 



(12) 



where Wij is a Gaussian having a 'width' of 2a — Is, 
truncated at ±3ct. 



(13) 



27rcr 



Simulated occultation lightcurves with no noise were 
normalized in the same way to evaluate the effects of 
smoothing on the events being sought. The amplitudes of 
these events were found to change by <0.5%, indicating 
that the normalization does not significantly attenuate 
or suppress an occultation signal. 

Figure |3] shows a segment of a time series for the 
eastern-most star in the M35 field. The smoothed time 
series is superimposed, and the normalized time series is 
shown in a separate panel. 
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Fig. 4. — A 10s sample of the time series for the eastern star in 
the M35 field. Both the raw flux counts (top) and the normalized 
data (bottom) are shown. The smoothed normalization curve is 
shown superimposed over the raw flux. Gaps mark the dead time 
of the camera. 



Imaging continued through periods of reduced trans- 
parency, or partial cloud. When a star's measured flux 
dropped by more than ~30% in a given time segment, 
the noise levels in that time segment increased accord- 
ingly by ~15% {Vn photon statistics), and the segment 
was masked. 
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The RMS noise values of the observed flux were com- 
puted for each target. The variates included all data for a 
given source; they represent the overall observation qual- 
ity. A 3(T-clipped standard deviation was used to avoid 
the influence of outliers. The eastern and western stars 
in the M35 field had noise levels of -4.3% and -3.4%, 
respectively. Both stars in the M34 field had noise levels 
of -1.6%. 

5. THE DETECTION ALGORITHM 

An object with a diameter >1 Fsu (1285m for visi- 
ble light at 40 AU) will produce a >40% flux decrease 
during occultation. No such events were found and a 
cross-correlation detection algorithm was used to search 
for events caused by sub-Fresnel-scale objects near the 
noise limit. 

Template diffraction patterns were cross-correlated 
with each time series to detect candidate occultation 
events. For a discretely-sampled time series, ti, and 
a template 'ke rnel', ki, the cross-correlation is (see 
lBracewelllll986l) : 

it*k),=J2uk,+,. (14) 

i 

A positive peak in (t ★ k)j indicates that t is similar to 
k when k is offset to position j. The detection threshold 
was set to > +6ati,k (random probability < 10^^). 

Masked data (gaps and high-noise regions) were ig- 
nored by the algorithm, and the total on-source observ- 
ing time is the number of non-masked points times the 
individual exposure time. 

To construct detection kernels for cross-correlation, 
theoretical occultation lightcurves were integrated and 
normalized in 0.025s segments to reproduce the effect of 
a camera taking exposures of that duration (our exposure 
time). Typically, 8-12 points were sufficient to model a 
kernel. 

At 40 Hz sampling, the diffraction features of a KBO 
occultation lightcurve are slightly under-sampled. A 
small shift in the sample start time produces subtle 
changes in the structure of the resulting kernel, and the 
cross-correlation could not be performed with kernels for 
all possible time-shifts. To construct detection-kernels, 
the integration start-time was synchronized with the 
event such that the kernel produced would be symmet- 
ric, and this single kernel was used as the the detection- 
kernel. Simulated occultation-kernels which were planted 
in a time series for testing purposes, were integrated with 
a random start-time. Simulated occultation events were 
'planted' in time series to evaluate the performance and 
recovery limits of the cross-correlation algorithm. For 
each planted event, a position in the time series was cho- 
sen randomly, and the relevant points in the time series 
were multiplied by the corresponding intensities from the 
kernel. 

5.1. Performance of the Detection Algorithm 

Artificial time series were used to assess the perfor- 
mance of the detection algorithm as they were certain 
to contain no occultations. Analysis of our time series 
showed that they were not composed of pure Poisson 
noise, but rather had power spectra with slopes of ap- 
proximately -1 (due to scintillation). Random variates 



were, therefore, assigned in the frequency domain to pro- 
duce power spectra with the appropriate slope, and these 
were inverse Fourier-tr ansformed to generate time series 
(jBickerton et al.ll2006D . This method of generating ar- 
tificial time series will be discussed in detail in future 
work. 

The signal-to-noise ratio (SNR) of the cross-correlation 
method was compared to a 'low-flux' method (a search 
based on identifying unexpectedly low flux values). We 
define 'noise' to be the root-mean-squared (RMS) devia- 
tion of the intensities from their mean, and 'signal' to be 
an event's maximum intensity deviation from the mean 
intensity. Planted occultation events for three different 
sized KBOs (at 40 AU) are shown in Figure[5]to illustrate 
the advantage of the cross-correlation method. 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 
Time (sec) Time (sec) Time (sec) 



Fig. 5. — The cross-correlation and low-flux detection meth- 
ods. Planted occultation events are shown for 4, 6, and Scr cross- 
correlation detections. Upper panels show the events in the time 
series data; the lower panels show the cross-correlation values. The 
left vertical axes are in units of standard deviation from the mean. 
The normalized flux units are shown on the right vertical axes in 
the upper panels. Original units for the cross-correlation are not 
shown. SNR values are higher for the cross-correlation variates 
than in the original time series. 

Any value above the detection threshold was compared 
to its expected value before being considered a candi- 
date. The expected value was estimated to be the max- 
imum value of the auto-correlation of the detection ker- 
nel ((fc -k k)j). High values in the cross-correlation were 
empirically found to be systematically smaller (0.93 x) 
than the auto-correlation value, and the auto-correlation 
value was scaled accordingly. To cover the range of possi- 
ble cross-correlation values, those within ±8a (more than 
sufficient to include all legitimate occultation candidates) 
of the scaled auto-correlation value were accepted as can- 
didates. 

The cross-correlation method yielded a higher signal- 
to-noise than the low-flux method for all objects tested 
(rA'BO = 150-500m, distance=40 AU). As even a conserva- 
tive (low) estimates of th e slope of the KBO size distribu- 
tion is quite steep (a=3. |Pan fc Sarill2005[ ). the ability to 
detect objects with radii only a few dozen meters smaller 
could significantly increase the total number of objects to 
which the observations are sensitive! 

The time series variates and those of the cross- 
correlation were consistent with having been drawn from 
Gaussian distributions (based on tests) . As we report 
no positive detections, the issue of false positive rates is 
not addressed here. It will be presented in future work. 
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5.2. Choosing the Detection Kernels 

The structure of an occultation lightcurve is influenced 
by the size of the KBO, the distance to it, and the min- 
imum distance between the KBO center and that of the 
stellar disk (called the impact parameter). Multiple de- 
tection kernels were used to provide coverage of the pa- 
rameter space. 

To select a kernel set, occultation events were planted 
in artificial data with a 4% noise level (the highest mean 
noise level in our M35 time series). With the distance 
and impact parameter held constant, objects of various 
sizes were planted and the detection algorithm was run. 
A detection kernel for a single object size was used to 
determine the range of object sizes it could successfully 
detect. This was repeated with detection kernels for 
different-sized objects, and a set with overlapping ranges 
of sensitivity was selected. This process was repeated to 
determine the ranges of sensitivity in terms of distance, 
and impact parameter. 

The final kernel set included objects at distances of 10, 
20, 40, 80, and 160 AU, and having sizes of rKBO=0-08, 
0.16, 0.24, 0.32, 0.40, and 0.48 Fsu (an event generated 
by an object larger than tkbo ^0.5 Fsu could be iden- 
tified by visual examination of the time series, even with 
noise levels >4%). Ten evenly-spaced impact parame- 
ters were included at each size/distance. Thus, a total of 
300 kernels were cross-correlated with each time series. 
The same kernel set was tested with 2% noise levels and 
performed equally well. 

5.3. The Recovery Limit and Apparent Shadow Size 

The smallest KBO radius to which the detection 
method is sensitive, rmin, and the maximum impact pa- 
rameter at which a KBO of this limiting size could be 
reliably detected, hmax, were determined by a plant and 
recover process. Occultation profiles were computed for 
the projected stellar radii and photometric passbands of 
the M34 and M35 targets, and were planted in artificial 
time series with noise levels of 2% (M34) and 4% (M35). 
The recovery fraction was determined for a range of KBO 
radii (20m intervals) and impact parameters (60m inter- 
vals) . 

For each radius/impact parameter pair (r,b), 100 
events were planted and the detection algorithm was run. 
A given (r,b) pair was considered detectable if its re- 
covery fraction was >80%, and the detectability limits 

ij^ rain -i^ max 

) were determined for occulting KBOs at dif- 
ferent distances (see Tabled]). 

6. RESULTS 

The algorithm flagged 177 events for review: 133 in 
the M35 data (27 per hour), and 44 in the M34 control 
data (21 per hour). To review each reported detection, 
the time segment surrounding it (±1.5s) was examined 
and compared with that of the adjacent star, and with 
the theoretical lightcurve of the detection-kernel. Raw, 
unreduced flux values of the two time series were also 
examined to provide an indication of any artifacts which 
may have been introduced by the reduction and normal- 
ization processes. 

The cross-correlation detection method and the 
subsequent review process found no viable occul- 
tation candidates in our time series. All flagged 
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The limiting TNO surface densities indicated by the null 
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photometric fluctuations were clearly present in both 
stars, and the events were easily identified as being spu- 
rious. Variations occurring on time-scales comparable to 
the duration of the smoothing kernel (1.0s) were only 
common in the more unstable M35 time series, and did 
not resemble occultations. Those occurring on shorter 
time-scales were common in all time series, and did in 
many cases, resemble occultations. Examples are shown 
in Figure m 



r=1200m d=160AU ip=760m r=300m d=40AU ip=1260m 




-1 -0.5 0.5 1 -1 -05 0.5 1 

Time (s) Time (s) 



Fig. 6. — Lightcurves of two false positive events from the M35 
data. The upper panels show the normalized flux for both stars 
(flux for the neighboring star is offset for clarity). The parameters 
of the kernel which detected the events are printed above each plot 
(r=KBO radius, d=distance, and ip=impact parameter), and the 
corresponding theoretical occultation profiles are shown above the 
stellar lightcurves. For each event, the raw flux values for the two 
stars are shown in the lower two panels. In each case, the coincident 
flux decrease observed in the neighboring star indicates that the 
event is inconsistent with interpretation as a KBO occultation. 

To evaluate the validity of the recovery tests, objects 
were added to the real, normalized, time series prior to 
running the detection algorithm. This was done to ensure 
that the manual review process did not alter the recovery 
rates significantly. The sizes and impact parameters of 
the events added and recovered are shown in Figure [71 
with the 80% contour from the 40AU theoretical recovery 
tests (Section I5.3P overlaid for reference. The theoreti- 
cal limit shown is based on 40AU tests, and the planted 
events shown are those with distances of 30AU < d < 
60AU. Some of the planted events are outside the 80% 



A search for sub-km KBOs 



7 



contour as a result of the broader distance range. 
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Fig. 7. — The recovered sizes and impact parameters of planted 
occultation events for 30-60AU objects. Squares represent planted 
occultation events; filled if recovered, and open if not recovered. 
The 80% cont our from the theoretical recovery tests for 40AU ob- 
jects (Section 15.31 1 is shown with a solid line. 



With this nuU result, we can place an upper limit on 
the surface density of TNOs. Taking a line of length 1 
arcsec perpendicular to the direction of motion, TNOs 
of density o arcsec"^, moving at v arcsec s~^, will cross 
the line at a rate av s~^. After t seconds, avt objects 
will have crossed the 1 arcsec line. For a star positioned 
somewhere on the line, the probability that x TNO line- 
crossing events will lie within 1 impact parameter, 6, of 
the star (along the line) is described by Poisson statistics 
with the rate 2bvta. 



P{x; b, V, t, a) 



{2bvt(TY 



-2bvta 



(15) 



As (To is assumed constant, only observations expected 
to have the same value for tJo (same ecliptic latitude) 
should be included. 

The product bvt represents an area swept through by 
the target star, and could be integrated over b{t)v{t)dt 
to account for variability in b and v, or to incorporate 
results from multiple stars. Our values of b and v were 
similar for both stars, and did not vary in time. The time 
t represented the sum of observing time on both stars. 

The 95% confidence limits for a occur when the prob- 
ability of detection is 0.05 (ie. the complement of equa- 
tionfTS]). With x — events, equationfTSlcan be reduced, 
and solved in terms of ctq: 



Co 



ln|l - P\ 
-2bvt 



(16) 



The limiting densities indicated by our null result are 
presented in Table [51 Values of r„ii„ (the size for which 
the result is relevant) and bmax were determined in the 
recovery tests described in Section 15.31 The relative ve- 
locities were computed as described in Section [2?2| and 
the time values are the total observed times (masked data 
not included). 

The limit for typical KBOs (40 AU, /3 « 0) was placed 
on the cumulative luminosity function (CLF) by convert- 
ing the object size (rmin) to an R magnitude (GOl) (as- 
sumes geometric albedo of 0.04, subscripts denote units). 



rriR- 25.9 + 2.5 log 



50 / 



100 



2r 



km 



(17) 



A radius of 430m corresponds an R magnitude of 
^35.3. Although the limit is placed on the cumulative 
luminosity function, equation 1161 assumes a uniform ob- 
ject size. This assumption is reasonable for even the 
most conservative (shallow) estimates of the slope of the 
size distribution, as they place the majority of detectable 
KBOs at the smallest size limit. 

7. EVALUATION OF REPORTED TNO OCCULTATIONS 

TNO occultations have be en rep orted by two separat e 
groups: iChang et all (|2006f) . and iRoaues et al.l (|2006f) . 
The TNO densities implied by these positive results are 
surprisingly high; here we evaluate their methods and 
results. 



7.1. TNOs detected bv \Chana etHU (MM) 

IChang et all (|2006[) (henceforth C06) searched -89 
hours of time-tagged x-ray photon data from the Rossi 
X-ray Timing Explorer (RXTE) archive, and reported 
58 TNO occultation events. Many of the events were 
coincident with high-energy particle arriv als recorded 
by ot her instrumentation on-board RXTE (jJones et al.l 
|2006| ). a nd 12 occultation c an didates remained a fter re- 
analysis [Gh^^eFaD (|2007l ). iJones et all (|2007| ) report 
that, at most, 10% of the observed lightcurve dips could 
be due to TNO occultations; and they suggest that it 
would be a mistake to conclude that any TNO occulta- 
tions have been detected in the Sco X-1 X-ray data. 

Our observations were not sensitive to objects in the 
size range observed by C06, and there is no disagree- 
ment with our null result. But the COG result indicates a 
surface density ~ 10^ times greater tha n that estimated 
by ext rapolation of the CLF reported bv lBernstein et aTl 
(^2001). 

The x-ray occultation candidates reported by C07^ 
were compared to model events (see Figure [S]). The Sco 
X-1 x- ray source was modele d with A=0.2-0.6nm light (6- 
2 keV. lBradshaw et al.ll2003[ ). produced by a circular disk 
with a diameter of 3m projected at 43 AU (C06). Occul- 
tation shadows were made for 25m, 50m, 75m, and 100m 
KBOs at 43 AU. The x-ray (A=0.4nm) Fresnel scale at 43 
AU is 35m and these shadows are expected to be domi- 
nated by diffraction effects. The 25m and 50m models are 
of order <1 Fsu and will produce circularly-symmetric 
diffraction shadows, regardless of the projected shape of 
the occulter. The 75m and 100m models correspond to 2, 
and 3 Fsu. Real (irregularly shaped) objects this large 
would produce a more complex pattern of fringes, but 
the overall width (duration) and depth of the lightcurves 
would be comparable to those modeled here with a cir- 
cular mask (the symmetry of the shadows, as it relates 
to the size of the occulter, is further discussed in Sec- 
tion[ni). A relative velocity of 25km s^^ (consistent with 
a Keplerian orbit at 43 AU) was used^^. 

The 12 reported events range in duration from 1.5ms 
to 2.0ms (3 or 4 low flux measurements), and exhibit 



^ Values for the flux levels were taken from lChane et aTl II2007I ). 
Figure 12 

^ Relative velocities were not provided by C07. 
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Fig. 8. — Reproductions of the 12 candidate occultation events 
(C07) compared to model lightcurves. The events virere manually 
reproduced from Figure 12 of C07, and event numbers are those 
assigned by C07. The model events simulate Sco X-1 occultations 
of 0.2-0.6nm x-rays by TNOs at 43 AU having diameters of: 25m, 
50m, 75m, and 100m. 

sharp transitions "when dropping from the baseUne in- 
tensity and "when returning to it. Sections through each 
of our modeled shadows were taken at an impact pa- 
rameter that maximized the depth-to-width ratio of the 
diffraction shadow (b=40m). Non-zero impact param- 
eters also avoided the flux increase associated with the 
central flash^ which was not observed in any of the C07 
events. 

Lightcurves for the 100m and 75m diameter TNOs are 
too broad (long in duration) to be consistent with the 
observed events, and the 25m object produces intensity 
changes which are too shallow. The 50m object's model 
lightcurve comes closest to reproducing the the observed 
flux changes. However, events with only 3 or 4 low flux 
measurements are too narrow to be consistent with oc- 
cultations by KBOs (C07 event numbers 4, 5, 14, 19, 52, 
59, 61, 66, 74, 90, and 104), particularly in cases where 3 
flux measurements are below ~40% of the original inten- 
sity (C07 events 4, 5, 52, and 66). The abrupt decrease 
in flux observed in the remaining event (event 41) and 
many others make them unconvincing when compared to 
the model lightcurves. 

The absence of the central flash in all events is statis- 
tically unusual. For objects of the purported size, the 
center of the shadow should display some rise in bright- 
ness, regardless of any irregularities in the shape of the 
occulter. 

As dozens of similar events were found to be instru- 
ment artifacts, acceptance of these detections as TNO 
occultations appears premature. 

7.2. Three TNOs detected bv lRoques etd\ (MM) 

iRoaues et ahl (|2006[ ) (henceforth R06) reported 3 TNO 
occultation events found in 604 minutes of 45 Hz imag- 
ing (simultaneous observations with Sloan g' and i' fil- 
ters). The reported objects had radii and distances of: 
110m at 15AU, 300m at 210AU, and 320m at 140AU, 
and would be beyond the detectable limits of our obser- 
vations. There is no disagreement with our null result. 

The sizes of the occulters present a serious inconsis- 

* For occulters smaller than 1 Fsu (or for a circular occulting 
mask of any size), the intensity at the center of the diffraction 
shadow is as though the occulter were not present. The feature, 
called a 'central flash', is due to constructive interference. 



tency when compared to the C06 result. With ~10 hours 
of data, sensitive to objects with impact parameters we 
estimate would be comparable to our own (b « 1500m, 
due to diffraction), R06 found two ~600m objects. The 
~89 hours of X-ray data from C07 would be sensitive 
to 600m objects with impact parameters of b «300m 
(no x-ray diffraction effects for an object this large) , and 
would be expected to yield ^ {bvt)cm/ Q}vt)m& = l-8x 
as many 600m objects - ^--^4. C07 reported no flux drop- 
outs consistent with objects this large. Though we found 
the x-ray flux drop-outs to be unconvincing as occulta- 
tions by 50m TNOs, events produced by 600m TNOs 
(had they been present) are unlikely to have been over- 
looked by C06/C07. Consider that at 140-210 AU, the 
x-ray (A=0.4nm) Fresnel scale is 65-79m - much smaller 
than the purported objects. Occultations would have 
produced near complete attenuation of the x-ray source 
(Sco X-1). At ~25 km/s relative velocity, a 600m object 
would occult the source for 24ms and would occupy 48 
consecutive photometric samples (C07 used 0.5ms sam- 
pling). Even an off-center chord through such a shadow 
would have produced an obvious occultation which would 
easily have been identified by the authors. The absence 
of any corroborating 600m TNOs in the x-ray time series 
is a serious inconsistency between R06 and C07. 

R06 use a 'variability index' (VI) for their detec- 
tion algorithm. The standard deviation was com- 
puted in 0.08-0. 40s intervals through the time se- 
ries (each denoted 'sig(int)'), and normalized by 
the number of standard deviations (stddev(sig)) from 
the overall mean standard deviation (meansig); thus: 
VI=(sig(int)-meansig)/stddev(sig). Localized increases 
above VI=4.5 were interpreted as being due to the bright- 
ness fluctuations associated with occultation candidates. 
Three events were identified with VI=5.3, 5.6, and 7.2. 
These are the moduH of the vector (Vlg/, Vli/), but as 
the Vli' values are small^, the modulus of the vector is 
well approximated by the VIg' value. 

The VI values were interpreted as (Vl)-sigma results, 
but the underlying distribution of standard deviations 
(sig(int) values) used to compute the VI values is non- 
Gaussian; it is positively skewed. 

To determine the VI values expected in the data, a 10 
hour (45 Hz) time series with 1.5% 1/f noise was artifi- 
cially generated (guaranteed to contain no occultations) . 
The VI was computed with intervals of 0.08s, 0.20s, and 
0.40, consistent with those described in R06. For each 
interval, the distribution of values larger than VI=4.5 
(the R06 threshold) are shown in Figure [51 The number 
of VI>4.5 values is higher than expected for Gaussian 
variates. Values as high as those observed by R06 are 
expected with this method. Values as high as VI=7.2 
were not observed, but a temporary increase in the noise 
level of a time series could easily produce such a mea- 
surement. 

In this evaluation, several tens of points were found 
with VI>4.5. This behavior is not reported by R06. 
However, the method produces a positively skewed distri- 
bution of variates which can be seen R06 Figure 1. The 
best that can be stated at this time is that the statistical 
significance of the points is lower than would be the case 

^ The events do not appear as strongly in the i' time series. This 
is expected as a result of the wavelength dependence of diffraction. 
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Fig. 9. — The distribution of VI values expected for 1/f noise. VI 
values were computed for intervals of 0.08s (bottom), 0.20s (mid- 
dle), and 0.40s (top); for a 10 hour (45 Hz) time series of 1.5% 
1/f noise. Dashed lines indicate the levels for Gaussian variates. 
High VI values occur more frequently than expected for Gaussian 
variates. 

for Gaussian variates, but is otherwise unclear. 

8. THE SURFACE DENSITY LIMIT 

The slope of the size-distribution for the sub-km ob- 
jects is estimated to be q~3.5. A variety of indepen- 
dent methods have been used, inclu ding: numerical sim- 
ulations of the coUision al evolution (|Kenvon &: BromlevI 
Pan & Sari'2005") , observations of cratering on Tri- 
ton (Stern fc McKinnon 2000), and inferences based on 
the absence of optical and infrared light pr oduced by such 
a population (jKenvon fc Windhorstll200ll ). 

Although we are skeptical of the reported sub-km TNO 
occultations for the reasons stated in Section [71 neither 
C06, nor R06, nor this work have reported any occulta- 
tion by a TNO larger than 1 km. At this limiting size, an 
upper limit more stringent than our own can be placed 
on the current sky surface density of TNOs. 

Equation \W\ was recalculated with the bvt products 
for each data set summed. Velocity and impact pa- 
rameter values of {v — 23km s^^ ,b — 1.5km) and 
{v = 25km ,b = 500m) were used for the R06 and 
COG data set, respectively. The velocities are the aver- 
age relative velocities for KBOs during the observation 
periods. The impact parameter for R06 (b= 1.5km) was 
chosen to accommodate diffraction effects. Such effects 
would be negligible for an x-ray occultation by a 1km 
object. 

Thus, the sum of bvt over each target is 
X^i ^i^i^i = 5.4 X 10~^° deg^, yielding an overall 
95% confidence limit on the surface density of 1 km 
(diameter) KBOs at 40 AU of ^ 7V(TOfl < 35.0) = 
2.8 X 10^ deg~^. This limit is shown in Figure [TOl 

9. OTHER CONSIDERATIONS 

Non-spherical objects originating in the outer solar sys- 
tem (eg. peanut-shaped comets) suggest that the spheri- 
cal symmetry assumed in our lightcurve models is not to 
be expected for TNOs. However, diffraction shadows cast 
by objects with Kkbo ^0.5 Fsu are circularly symmet- 
ric, regardless of the shape of the occultcr. Occultation 
shadows for different-sized elliptical TNOs {a/b = 2.5, 
where a/b is the ratio of the semi-major and semi-minor 
axes) were generated to illustrate this (see Figure [TT|). 
The circular symmetry of the shadows is preserved for 
aKBO <0.5 Fsu (-650m for visible light at 40 AU), de- 
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Fig. 10. — The CLE showing the overall upper limit on 1km 
KBOs. Upper limits for 95% and 50% confidence are shown with 
horizontal lines above a downward pointing arrow. Slopes of a=0.6, 
0.7, a nd 0.8 are shown with the model CLF of Bernstei n et al.l 
II2004I ) . A vertical line is used to indic ate the current limit for 
direct observation IIBernstein et al]|2004l) . 

spite the asymmetry of the occulting masks. The shadow 
profiles for objects larger than this do not require cross- 
correlation to detect; they would have been identifiable 
by visual inspection of the time series, had any been 
present. 
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Fig. 11. — Occulting masks (upper panels) and diffraction shad- 
ows (lower panels) for different-sized non-spherical KBOs (a/b=2.5 
in each case, where a and b are the KBO's semi-major and semi- 
minor axes.) The diffraction shadows are circularly symmetric for 
a<0.5 Fsu (~650m for visible light at 40 AU), despite the asymme- 
try of the occulting masks. The grayscales were chosen individually 
to cover the range of each shadow's intensity. 

Plant and recover tests were performed with profiles 
for 500 elliptical objects having uniformly distributed 
random parameters of: 200m< a < 1000m, Om<impact- 
parameter<4000m, 1.0< (a/b) <3.0, and < 9 < n; 
where a, b, and 9 are the semi-major and semi-minor 
axes of the object, and 9 is the position angle of the 
rotated object. Kkbo was taken to be 0.5(a -I- b) in or- 
der to compare the results to Figure [T] The observed 
detection limits did not change and are shown in Fig- 
ure 1121 Thus, cross-correlation detection kernels based 
on circularly-symmetric occulting masks successfully de- 
tect occultations by sub-Fresnel-scale irregularly-shaped 
KBOs. 

Our method of generating artificial time series did not 
incorporate the effects of variable noise levels in a time 
series. Although many of our recovery tests were per- 
formed with constant-noise artificial time series, they 
were verified in our final detection run by planting artifi- 
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Fig. 12. — The recovered sizes and impact parameters of planted 
occultation events for elliptical objects. Squares represent planted 
occultation events; filled if recovered, and open if not recovered. 
The 80% contour from th e th eoretical recovery tests for circu- 
lar 40AU objects (Section I5.3|l is shown with a solid line. The 
planted events had parameters of: 200m< a <1000m, Om<impact- 
parameter<4000m, 1.0 < (a/b) < 3.0, and < 8 < it; where a, b, 
and 6 are the semi-major and semi-minor axes of the object, and 
6 is the position angle of the rotated object. Kkbo was taken to 
be 0.5(a + fe). 

cial occultation events in the real data. The results were 
found to be consistent with simulations we performed 
with the artificial data, and we believe that the presented 
estimates for the limiting detectable KBO size, rmm, and 
the maximum impact parameter, bmax, are accurate. 

The presence of a period of increased noise within a 
time series alters the statistical properties of the variates. 
The overall measured standard deviation of the variates 
will lie between those measured for the low-noise and 
high-noise segments, which will lead to increases in the 
rates of false negatives and false positives in the respec- 
tive segments. Our time series were broken into segments 
having approximately constant noise levels to mitigate 
this problem. 

10. SUMMARY 

Analysis of our 5.0 star-hours of photometric time se- 
ries on two B9V stars in the ecliptic open cluster, M35, 
revealed no viable candidate occultations. For the first 
time, the recovery of artificial occultation lightcurves 
was used to determine the limiting detectable object 
size and impact parameter for a KBO occultation. We 
place an upper limit on the surface density of KBOs of 
J^NirriR < 35.3) = 3.5 x 10^° deg-^ (m^ -35.3 corre- 
sponds to a ~860m (diameter) object with a geometric 
albedo of 0.04). 



Having e valuated positive detections of TNO occulta - 
tion events (|Chang et al.l[2006l . [20071 : [Rogues et al.l2006h . 

we believe that the published events are likely to be spu- 
rious and that no serendipitous stellar occultation by a 
TNO has yet been reported. Comparison of the expected 
lightcurve structure for the C07 events with the observed 
structures show that the decreases in flux are too brief 
to be consistent with occultations by small bodies at 40 
AU; and evaluation of the variability index, VI, used by 
R06 revealed that the VI variates are not normally dis- 
tributed, but are drawn from a positively skewed dis- 
tribution, prone to yielding high values. Also, the two 
~600m (diameter) TNOs discovered in 10 star-hours of 
R06 data suggest ~4 such events should be present in the 
89 star-hour x-ray time series observed at a similar eclip- 
tic latitude i(3Roe = OA°S, OA°S, 7.1°iV, /3co6 = 6°iV). 
Nothing this large was reported by C06/C07. 

Combining all results from the different surveys we 
place an overall 95%-confidence upper limit on the TNO 
surface density of: J^^irnR < 35.0) = 2.8 x 10^ deg-^. 
A direct extrapolation of the observed CLF slope 
{a wO.7) would yield a density of J^^i^R < 35.0) = 
1.9 X 10^ deg-2. 

Our work is the first complete and careful examina- 
tion of observational data for KBO occultations and their 
proper analysis and interpretation. 
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